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Abstract 

We first discuss the geometrical construction and the main mathe- 
matical features of the maximum-entropy-production/steepest-entropy- 
ascent nonlinear evolution equation proposed long ago by this author in 
the framework of a fully quantum theory of irreversibility and thermo- 
dynamics for a single isolated or adiabatic particle, qubit, or qudit, and 
recently rediscovered by other authors. The nonlinear equation generates 
a dynamical group, not just a semigroup, providing a deterministic de- 
scription of irreversible conservative relaxation towards equilibrium from 
any non-equilibrium density operator. It satisfies a very restrictive stabil- 
ity requirement equivalent to the Hatsopoulos-Keenan statement of the 
second law of thermodynamics. We then examine the form of the evolu- 
tion equation we proposed to describe multipartite isolated or adiabatic 
systems. This hinges on novel nonlinear projections defining local opera- 
tors that we interpret as "local perceptions" of the overall system's energy 
and entropy. Each component particle contributes an independent local 
tendency along the direction of steepest increase of the locally perceived 
entropy at constant locally perceived energy. It conserves both the locally- 
perceived energies and the overall energy, and meets strong separability 
and non-signaling conditions, even though the local evolutions are not 
independent of existing correlations. We finally show how the geomet- 
rical construction can readily lead to other thermodynamically relevant 
models, such as of the nonunitary isoentropic evolution needed for full 
extraction of a system's adiabatic availability. 

Keywords: entropy, irreversibility, nonlinear quantum dynamics, steepest en- 
tropy ascent, maximum entropy production principle, quantum thermodynam- 
ics, Onsager reciprocal relations. 
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1 Introduction 



In this paper, we construct a class of model evolution equations (applicable not 
only to open systems but also to closed isolated systems) capable of describing — 
simultaneously and in competition with the usual Hamiltonian unitary evolution — 
the natural tendency of any initial nonequilibrium state to relax towards canoni- 
cal or partially-canonical thermodynamic equilibrium, i.e., capable of describing 
the irreversible tendency to evolve towards the highest entropy state compatible 
with the instantaneous mean values of the energy (and possibly other constants 
of the motion and other constraints). 

In this introduction, we briefly review some essential challenges of a prevail- 
ing model of irreversibility (Section ll.ip . we briefly discuss the original motiva- 
tion that lead us to develop a quantum maximal entropy production formalism 
(Section 1 1 - 2[) . we discuss the geometrical framework of our derivation (Section 
II. 3[) . and we outline the structure of the rest of the paper fSection ll.4[) . 

1.1 Kossakowski-Sudarshan-Gorini-Lindblad quantum mas- 
ter equation 

The prevailing model of irreversibility starts from unitary dynamics but assumes 
that no system is truly isolated, so that even an initial pure state becomes mixed 
due to increasing system-environment entanglement. 

The system-environment entanglement builds up due to interactions accord- 
ing to the standard Liouville-von Neumann unitary dynamics of the overall 
system-environment composite. In this phenomenological model a system A is 
assumed to be weakly coupled with a reservoir R, so that they can exchange 
energy via unitary evolution of the overall state par- The reservoir R is mod- 
eled as a collection of a large number of quantum systems (many degrees of 
freedom, e.g., the modes of the electromagnetic field). Because of the weak 
coupling, the unitary dynamics of par produces both an energy exchange and a 
build up of correlations between the system and the reservoir. However, justified 
only by some heuristic reasoning, a crucial additional assumption is injected in 
the derivation (Markovian approximation): that correlations smear out rapidly 
enough so as to maintain A and R effectively decorrelated not only initially, 
but at all times. One rationale usually offered for this approximation is that 
when the reduced density operator of A is time averaged (coarse grained) over 
a sufficiently long time interval, which is nevertheless still much shorter than 
the system's time scale of interest, the average correlations becomes negligible, 
and the averaged state effectively factors at all times. The model is phenomeno- 
logical and basically charges the reservoir's complexity for the system's losing 
quickly its memory of past interactions^ By tracing out all the reservoir's de- 

The literature on the well known problem of accounting for the arrow of time and the origin 
of irreversibility within Statistical Mechanics and its time reversible underlying dynamics is 
too vast to adequately review, see, e.g., [I]. Many alternatives to time averaging have been 
proposed as rationales for coarse graining (see, e.g., [2] [3j H] ) • We agree with a Referee, that 
a common feature of all these attempts to justify the continuing validity of the KSGL master 
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grees of freedom, the overall unitary dynamics gives rise (under such Markovian 
approximation) to a system's reduced dynamics which is nonunitary, linear, 
completely positive and generated by the celebrated Kossakowski-Sudarshan- 
Gorini-Lindblad (KSGL) quantum master equation [Sj 

| = -jr[B*P] 4E( 2 ^ - {V}Y,,p}) , (1) 

3 

where the Vj's are some operators on 7i (each term within the summation, 
often written in the alternative form [Vj, pVj] + [Vjp, Vj], is obviously traceless). 
It has been used for a number of successful models of dissipative quantum 
dynamics of open subsystems. Operators Vj in |T]) are in general interpreted 
as creation and annihilation, or transition operators. For example, by choosing 
[6], Vj = Cr S \r)(s\, where c rs are complex scalars and |s) eigenvectors of the 
Hamiltonian operator H, and defining the transition probabilities w rs = c rs c* s , 
the KSGL equation becomes 

^ = -j[H,p} + J2 w rs(\s)(r\p\r){s\- ±{\ S )(s\,p}) , (2) 



dt h L ' rj ^— ' y 2 

rs x 

or, equivalently, for the nm-th matrix clement of p in the H representation, 
dp m 



df 



which, for the n-th energy level occupation probability p n — p nn , is the cele- 
brated Pauli master equation 

~~j~T = w nrPr ~ Pn ^ Wrn ■ ( 4 ) 

r r 

In this widely accepted model, the assumption of erasure of correlations is 
the sole mechanism responsible for "entropy generation" [7], but the apparent 
contradiction with the assumed underlying reversible unitary dynamics, i.e., 
the Loschmidt paradox is still lurking behind. The situation is parallel to what 
is needed to "derive" the classical Boltzmann equation from the underlying 
reversible Hamilton-Liouville dynamics. 

Philosophically, we find it hard to understand how diffusion of mass, momen- 
tum, energy, and charge, could find their justification in a "loss of information 
on the time scale of the observer leading to rapid decoherence from the en- 
tanglement which continuously builds up by weak coupling with environmental 
degrees of freedom". Is this the real physical reason for the "universal ten- 
dency in nature to the dissipation of mechanical energy" already recognized by 
Kelvin in 1852 [5]? Do we have alternatives to understand and model physical 
phenomena that are manifestly time asymmetric? 

equation, or equivalents, is that they require additional, question begging assumptions that 
are at odds with an overall unitary evolution. 
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1.2 Locally maximal entropy production dynamics as a 
nonlinear alternative to the KSGL model equations 

With these kind of motivation, thirty years ago we designed a possible alter- 
native based on the assumption, still to be validated or invalidated, that irre- 
versibility is a fundamental microscopic dynamical feature and as such it must 
be built into the fundamental laws of time evolution. Therefore, we constructed 
a fundamental non-unitary extension of standard Schrbdinger unitary dynamics 
not contradicting any of the successful results of pure-state quantum mechanics, 
and yet entailing the Second Law as well as an objective entropy increase for 
mixed states of an isolated system. We have shown in Refs. [SJ QH [HI HI] 
that such an approach is possible based on a steepest entropy ascent, i.e., max- 
imal entropy generation, nonlinear and non-unitary equation of motion which 
reduces to the Schrodinger equation for pure states. A 1985 Nature editorial [13] 
defined this approach "an adventurous scheme which may end arguments about 
the arrow of time" . Until now, however, the theory has been rather ignored 
and neither validated nor invalidated experimentally. Therefore ( "the proof of 
the pudding is in the eating" [13]), it remains just an interesting but little ac- 
knowledged and pursued theoretical alternative to the standard model. Some 
recent discussion about it, is found in [14] . The challenge with this approach 
is to ascertain if the intrinsic irreversibility it implies at the single particle (lo- 
cal, microscopic) level is experimentally verifiable, or else its mathematics must 
only be considered yet another phenomenological tool, at the same level as the 
quantum Markovian master equations which, as we have seen, are not free of 
their own challenges. 

The central conceptual difference between the proposed approach, and the 
approaches based on attempting to derive the KSGL equation, is that this ap- 
proaclQ regards a non-pure density matrix as representing a real ontological 
object, the actual state of the world, and is not understood as just an epistemic 
ignorance of which particular pure state the world is 'really' in|f] 

Therefore, we avoid the (unnecessary) severe restrictions imposed by lin- 
earity on the evolution equation, and we open up our attention to evolution 
equations nonlinear in the density operator p. It may at first appear natural to 

2 To our knowledge, the first pioneering work where this assumption is made explicit and 
used consistently to build a unified quantum theory of mechanics and thermodynamics, is 
Ref . |15| . We thank a Referee for suggesting the wording of this paragraph and the following 
footnote. 

3 If one assumes that the 'true' state of the world is actually a pure state, and that the 
probabilities in a density matrix can only be a reflection of uncertainty, then it is fairly easy 
to see that the evolution of the density matrix must be linear: if pi — > and p2 —* p' 2 , 
then a probabilistic mixture pipi + P2P2 — * Plp'i + P2P2- This linearity does not need to 
hold for our proposed evolution law as a 'real' density matrix p is a distinct physical state, 
even though numerically it may be equal to p = wipi + W2P2 with w\ = pi and U12 = P2, 
and is quite different from a probabilistic mixture arising through ignorance or uncertainty. 
Even if pi — » p' x and p2 —* p 2 > does not follow that p = wipi + W2P2 —* «J1Pi + W2P 2 - 
This conceptual difference is at the heart of our original approach. The physical reality 
attributed to the density matrix also legitimates treating the entropy — fc B Tr(plnp) as a 'real' 
physical quantity, in the manner of energy or mass, and not as an expression of information 
or uncertainty about a probability distribution. 
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maintain the Kossakowski-Lindblad form (fTJ) and simply assume that operators 
Vj are functions of p. This is true only in part for the evolution equation we 
will construct. Indeed, our hermitian operator AM in our Eq. below, can 
always be written as — ^ Vj(p)Vj(p) and therefore our anticommutator term 
may be viewed as a generalization of the corresponding term in ([1]) . 

However, our geometrically motivated construction based on the square- 
root of the density operator effectively suppresses the term corresponding to 
J2j Vj pVj m ©■ The reason we find this suppression desirable is the following. 

Due to the terms vj pVj ([I}, whenever p is singular, its zero eigenvalues may 
change at a finite rate. This can be seen clearly from Qj, by which dp n /dt is 
finite whenever there is a nonzero transition probability w nr from some other 
populated level (p r ^ 0), regardless of whether p n is zero or not. When this 
occurs, for one instant in time the rate of entropy change is infinite, as seen 
clearly from the expression of the rate of entropy change implied by ([T]) , 

^ = fc B ]T Tr(VjV jP ]np- VjpVj bxp) = k B ^(^rO^nrCPr - Pn) hxp r , 
j jrn 

(5) 

where p r denotes the r-th eigenvalue of p and (Vj) nr the matrix elements of Vj 
in the p representation. 

We may argue that an infinite rate of entropy change can be tolerated, 
because it would last only for one instant in time. But the fact that zero eigen- 
values of p in general would not remain zero (or close to zero) for longer than 
one instant in time, to us is an unphysical feature, at least because it is in 
contrast with a wealth of successful models of physical systems in which great 
simplification is achieved by limiting our attention to a restricted subset of rel- 
evant eigenstates (forming a subspace of Tt that we call the effective Hilbert 
space of the system [16]). Such common practice models yield extremely good 
results, which, being reproducible, ought to be relatively robust with respect 
to inclusion in the model of other less relevant eigenstates. In fact, such added 
eigenstates, when initially unpopulated, are irrelevant if they remain unpopu- 
lated (or very little populated) for long times, so that neglecting their existence 
should introduce very little error. The terms VjpVj, instead, would rapidly 
populate such irrelevant unpopulated eigenstates and void the validity of our so 
successful simple models. Of course, we may deliberately overlook this instabil- 
ity problem by making highly ad-hoc assumptions, e.g., by forcing the Vj's to 
be such that (Vj) nr = whenever either p n = or p r = 0. But, in this case, we 
can no longer claim true linearity with respect to p. 

Another important general physical reason why we find it advantageous that 
our construction excludes KSGL terms that generate nonzero rates of change 
of the zero eigenvalues of p, is that such terms are construed so as to pre- 
serve the positivity of p in forward time, but in general they do not maintain 
it in backward time. Such mathematical irreversibility of the Cauchy problem 
is often accepted, presented, and justified as a natural counterpart of physical 
irreversibility. However, we already noted in [16 that it is more related to a 
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principle of causality than to physical irreversibility. The strongest form of a 
general non-relativistic principle of causality — a keystone of traditional physical 
thought — requires that future states of a system should unfold deterministically 
from initial states along smooth unique trajectories in the state domain, defined 
for all times, future as well as past. Accepting mathematical irreversibility of 
the model dynamics, would imply giving up such causality requirement. The 
foundational virtue of our dynamical group is in its very existence, which shows 
a simple conceivable alternative whereby we are not compelled to cope with 
such a major conceptual loss. Regardless of these important but highly contro- 
versial foundational implications, we have shown in [17] that our Eq. (|32[) can 
effectively describe relaxation within an isolated system, and yet it is mathe- 
matically reversible, in the sense that it features existence and uniqueness of 
well-defined solutions both in forward and backward time. 

Eq. (f32|) describes physically irreversible time evolutions, in the sense that 
the physical property described by the entropy functional — k B Tr(phip) is a 
strictly increasing function of time for all states except the very restricted set 
of equilibrium states and limit cycles defined by Eq. ([i3"j) below. 

Similarly to our presentation of MEPP dynamics in a general probabilistic 
but non-quantal framework in [181 119] , in this paper we focus on the mathe- 
matical features and the potential phenomenological applications of a quantal 
MEPP dynamical equation. We emphasize that the formalism has an intrin- 
sic mathematical validity per se as a mere phenomenological tool. This may 
be useful also for those who remain understandably skeptical about the cited 
adventurous scheme of our original attempt to unify mechanics and thermody- 
namics, whereby, again, we proposed a resolution of the long-standing dilemma 
about the arrow of time based on building the Hatsopoulos-Keenan statement 
of the second law directly into the dynamical postulate of quantum theory. 

On the other hand, if proved valid at the fundamental level as envisioned 
in its original framework, our nonlinear dynamical law would imply the incom- 
pleteness of unitary pure-state zero-entropy quantum mechanics and the need to 
broaden it as suggested in [151 [20j . In such context, our microscopic dynamical 
theory might also be seen to accomplish the program sought for with limited 
success in the 1980's by the Prigogine school [HI [55], namely, to build a math- 
ematical theory of microscopic irreversibility (the question: "minimal entropy 
production or maximal entropy production?" is an interesting one, and, at least 
mathematically, is clarified in [23[ 124j where it is shown that maximal entropy 
production in general, implies minimal entropy production at some constrained 
stationary states). 

1.3 Maximal Entropy Production Path in a Maximum En- 
tropy Landscape 

The determination of a density operator of maximum entropy subject to a set of 
linear constraints has applications in many areas of quantum physics, chemistry, 
information, and probability theories 25, 26, 27, 28, 29 . The maximum entropy 
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density operator typically represents a thermodynamic equilibrium state or a 
constrained-equilibrium state of the system under study. 

Having set aside the cited implications on conceptual and physical quantum 
foundations, this paper focuses on the geometrical construction of our MEPP 
nonlinear quantum master equation, presented as the mathematical generaliza- 
tion of the maximum entropy problem to the nonequilibrium domain, by dis- 
cussing a general rate equation for the description of smooth constrained relax- 
ation of arbitrary non-equilibrium density operators towards maximum entropy. 
The nonlinear rate equation keeps the constraints constant at their initial values 
and increases the entropy until an unstable or stable maximum-entropy equilib- 
rium state is approached. The unstable equilibrium density operators are those 
with at least one zero eigenvalue and all others canonically distributed (see Eq. 
(|43|) below). The rate equation is also consistent with an Onsager reciprocity 
theorem interestingly extended to the entire non-equilibrium domain. 

Geometrically, it has a clear representation in square-root density operator 
state space. Every trajectory unfolds along a path of steepest entropy ascent 
compatible with the constraints (constrained geodesies). For an isolated system, 
the constraints represent constants of the motion. For more general quantum 
thermodynamics modeling, such as for rate-controlled constrained equilibrium 
modeling of chemical kinetics [551 US] , the constraints may be assigned a specified 
time-dependence. 

The well-known maximum entropy problem which sets our context (land- 
scape, to use the terminology of nonlinear optimization) is that of seeking a 
density operator p whose entropy S(p) = — fc B Trplnp is maximal subject to 
given magnitudes (Ak) of one or more linear constraints TrpAk = (Ak) for 
k = 0,1,..., n where Ak is the hermitian operator associated with the fc-th 
constrained observable. We assume the first constraint to be the normalization 
condition, so that Aq = I and (Ao) = 1. Moreover, as suitable to model a 
canonical isolated system, below we will assume for simplicity a single nontriv- 
ial constrained observable, the energy, represented by the Hamiltonian operator 
A X = H, 

The maximizing density operator p* can be written as p* — exp (— 53fc=i ^fc ^fc) 
with Q — Trexp(— X)fe=i ^kAk) where the Lagrange multipliers A& are deter- 
mined by the values (Ak) of the constraints. 

In this landscape, we wish to consider the following general problem in the 
non-equilibrium domain. We seek a time-dependent density operator, namely, 
an operator function (one-parameter family) p(t) , whose zero eigenvalues remain 
zero at all times and whose entropy S (p(t)) is maximally increasing with time 
t, i.e., 

max — k B — Ttp(t) In p(t) > subject to Tr p(t)Ak — const and ^ = const, 
at at 

(6) 

where £(t) is a properly formulated measure of length of a trajectory in density 
operator state space. In time interval dt, among all the possible trajectories in 
state space that have length d£, the system selects that which yields the maximal 
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increase in the value of the entropy functional, i.e., the path of steepest entropy 
ascent. This models a most irreversible quantum evolution towards maximum 
entropy. It is a realization at the (fundamental?) quantum dynamical level of 
the "principle of maximal entropy generation" [301 131] . The empirical validity 
of such a principle at the phenomenological level has been recently affirmed 
(explicitly or implicitly) by various authors in different fields and frameworks 
(see, e.g., [23l [32l [33l EH [33 [36] ) . 

Though overlooked even in recent reviews [36] . our nonlinear dynamics is 
one of the earliest instances and implementations of the ansatz which today 
goes under the name of maximum entropy production principle (MEPP). A 
reason for the oversight may have been that starting in 1984 (3U1 GH] we classi- 
fied our approach as 'steepest entropy ascent' dynamics, rather than MEPP, to 
emphasize that the qualifying and unifying feature of this dynamical principle 
is the direction of maximal entropy increase rather than the rate at which a 
noncquilibrium state is attracted in such direction. 

The formalism presented here has mathematical features of great generality, 
and is presented in a form readily adaptable to different applications. It was 
originally "designed" by the author in 1981 [6] and subsequently developed 
[HI [13 [HI HZl GSl IE] to obtain an equation of motion for a quantum theoretical 
unification of mechanics and thermodynamics [131 Q~5] . 

Recently, the original equation has been partially rediscovered in the same 
context [37) . The idea of steepest-entropy-ascent time evolution of a probability 
distribution has also been recently rediscovered in [38] but with important differ- 
ences we discuss in [17! . Because of its intriguing general features, we suggested 
long ago [18j that the formalism maintains its appeal even when abstracted from 
its original physics purpose because it provides a powerful mathematical tool 
for phenomenological modeling applications. It has been recently rediscovered 
also in such a broader maximum entropy formalism, probabilistic context |39j . 

1.4 Outline of the Paper 

In Section [2] we present the geometrical reasoning that leads to the construction 
of our main equation, Eq. 1321 for a single particle system. Important to this 
development (as well as the rest of the paper) is the material in the Appendix, 
which reviews some well known but little used geometrical notions, and sets the 
notation of our derivations. 

In Section [3] we outline the main features and theorems of Eq. [32] 
In Section [4] we discuss the generalization of our dynamics to composite 
systems, which is nontrivial in view of the nonlinearity of the steepest entropy 
approach. 

Finally, in Section [5] we discuss a further generalization of the foregoing non- 
cquilibrium problem whereby the magnitudes (Ak) of the constraints and the en- 
tropy rate of change may be assigned definite or interrelated timc-dcpcndcnces. 
This may become useful in the framework of quantum thermodynamics mod- 
eling of a non-work interaction, by which we mean [40 an interaction where 
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in addition to energy exchange between the interacting systems, there is also 
entropy exchange. 



2.1 



Geometrical construction of a single-particle 
MEPP quantum dynamics 

Reformulation in terms of square-root density opera- 
tors 



Because we seek a well-defined time evolution equation for the density operator, 
we must enforce at all times the positive semi-definiteness and hermiticity con- 
straints, i.e., p(ty — p{t) > 0. To this end it is convenient to change variables 
and represent quantum states by means of the square-root density operator 7 
defined as follows 



7 = EVP. 



w = u~ 



7 f 7 



(7) 



where ^fp is the positive square root of p, and U an arbitrary unitary operator 
that in the end will turn out to be irrelevant, much like phase factors in usual 
quantum mechanics. In the original derivations we assumed U = I, but as 
suggested in [37] the introduction of U has some formal advantage. 

Notice that were it only to cope with the positive semi-definite constraint, 
we could choose as new 'variable' any function of p whose inverse is even. The 
main reason for choosing ^fp is geometrical and part of the steepest entropy 
ascent assumption. We note here that, like done in Ref. [38] in a non quantal 
context, we could derive a steepest entropy ascent dynamics without switching 
to a square root representation, but in such case the entropy gradient would 
not be well-defined on the entire domain, it would diverge whenever one of the 
eigenvalues of p is zero, and as a result the dynamics would exhibit unphysical 
infinite-rate effects and would not conserve positivity when solved backwards in 
time. 

In order to introduce the geometrical notion of steepest entropy ascent, we 
need to define what we mean by "distance" between two density operators and 
by "length" of a portion of trajectory in state space, i.e., a one parameter family 
of density operators, a time evolution. The proper unique natural metric for this 
purpose is known in statistics as the Fisher-Rao metric (see e.g. [4T]l42. 43 ). For 
a one-parameter family of discrete distributions, p(t), where t is the parameter, 
the distance between distributions p(t + dt) and p(t) is 



\ 



dt 



x Y,-(^r) dt = 



\ 



j2(^E) ,/, 



dt 



= (i*f dt = Vx-xdt . (8) 

Thus, square-root probabilities Xi — ^ppl are the most natural variables in that: 
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• the space becomes the unit sphere, x • x = 1 (J2iPi = l)j 

• the Fisher- Rao metric simplifies to dl = Vx • x dt, or equivalently d£ 2 
tlx ■ cfoc; 

• the distance between any two distributions is the angle d(xj , x 2 ) = cos -1 (xi • 
x 2 ). 

We therefore conveniently rewrite the density operator formalism in terms 
of the square-root-density operator representation of states. To do so, we equip 
the space of linear (not necessarily hcrmitian) operators on 7i with the real 
scalar produclQ 

X Y = %Tr(X^Y + Y j X) . (9) 

The state space becomes the unit sphere 7 • 7 = 1 (= Trp with p = 7^7 au- 
tomatically positive semidefinite). On the state space we therefore adopt the 
Fisher-Rao type of metric ^(71,72) = cos~ 1 (7i • 72), so that along a time de- 
pendent trajectory 

d£ = 2y^r^dt. (10) 

2.2 Notation. Gradients of the Energy and Entropy Func- 
tional 

For simplicity, in addition to the normalization constraint, here we will assume 
a single additional constraint, namely, energy conservation^ with associated 
Hamiltonian H hcrmitian on TL. The extension to more constraints is straight- 
forward [44] in view of the formalism in the Appendix. In terms of the square- 
root density operator 7, the functionals representing the mean values, their time 
rates of change along a time dependent trajectory j(t), and the dispersions and 
covariance of the energy and the entropy, are conveniently rewritten introducing 

4 Note that this real inner product on the vector space of linear operators does satisfy the 
necessary rules, including of course that X-Y = Y- X,X-X>0, and X ■ X = iff X = 0. 
It clearly differs from the more usual complex inner product Tr(Xty). 

5 If H depends on aset of time dependent parameters Ai, A2,. . . , \k > the system is adiabatic 
and we usually interpret (dH/dt) = ~^2j,(dH/d\k)^k as the rate of work exchange between 
the system and a set of K "work elements" 1151 mechanically coupled with the system through 
the variation of these parameters. Such an adiabatic system undergoes what in 1401 we call 
a "weight process", and the energy balance equation (energy conservation) reads d{H)/dt = 
(dH/dt) or, in the notation we introduce in this section, H' '7 = 0. 
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the following notation^] 



±1 


— A j 11 , so inai 11 yll — \11 ) — 2 M 




o 


— 9^/<? with Pr> Inn — — h n Pv \ m 

— ^("i Willi O — foB-i Ranp m p — K B Hl^/J -|- i^Kcip), SO 


that ("19"! 
nidi 1 1 




= — K B Lrpmp = 2*5 -7. 


/ 1 Q\ 
(13) 


(AH)' 


= 27A.fr, with AH = H - (H)I, so that (AH)' -7 = 0; 


(14) 


(AS)' 


= 2 7 AS*, with AS = S - (S)I, so that (AS")' -7 = 0; 


(15) 


P 


= 7^7 + 7^7, since p — 7^7; 


(16) 


Trp 


= 7'7 + 7- 7 = 27-7 = 0, since Trp = 7 • 7 = 1 ; 


(17) 


d(H)/dt 


- (dH/dt) = TrpH = 4 (7 • H' + 7 • H 1 ) = H' ■ 7; 


(18) 


d(S)/dt 


= /s B Trplnp)/di = £' • 7; 


(19) 


(AHAH) 


= Tip(AH) 2 = 7 Ai? • 7 Ai? = J (AH)' ■ (AH)'; 


(20) 


(ASAS) 


= Trp(AS*) 2 = 7 AS* • 7AS = I (AS)' ■ (AS)'; 


(21) 


(ASAH) 


= (AH AS) = i Tip{AS, AH} = jAS ■ -/AH = J (AS)' 


• (A# )'; 






(22) 



where Pa an p = 7 — -Fkerp is the projector onto the range of p (the subspace given 
by the eigenvectors of p with non-zero eigenvalues). Because of Eqs. (fTT|) , (fT%j) 
and (fi"9|) we call the operators 27, 77' and S' the normalization gradient, energy 
gradient and entropy gradient operators, respectively. In the same sense, (AH)' 
and (AS 1 )' are the gradients of the null functionals Tip(AH)' and Trp(A5)'. 

It is noteworthy that a dual set of definitions may be constructed using 
7 = JpU = 7 t, p = r/t, H' = 2^H, d(H)/dt - (dH/dt) = H' ■ 7 t, and so on, 
leading however to the same results. 

6 It hinges on the inner product denned by The logic and some details are as follows. 
For any hermitian A, we define A' = 27 A so that (A) = Tr(pA) = J A' • 7 (clearly, in 
general A' is not hermitian). When the ' operation is applied on A A = A — (A) I, we 
obtain (AA)' = 27AA and, in general, (AA)' -7 = 0. Next, because p = 7T7 and Trp = 
7-7 = 1, we have Trp = 7 ■ 7 + 7 ■ 7 = 27-7 = 0, 7 A ■ 7 = 7 A ■ 7 and, therefore, 
Tr(pA) = Tr[(7t 7 + 7+7)/!] = JTrftt^A) + (27A) t -y] = A' ■ 7 = (AA)' ■ 7. In general, 
therefore, d(A)/dt — (dA/dt) = Tr(pA) = A' • 7. For a time independent Hamiltonian H, 
d(H)/dt = H' ■ 7. Moreover, because Tr(pS) = (proof in the next footnote), we have 
d(S)/dt = S' • 7 in spite of S = — feB^Ranplnp being time dependent. Finally, for any 
hermitian pair A and B, we have the identity 7AA ■ 7AP = J ((AAAB) + (ABAA)). Because 
(AAAP) - (APAA) = Tr(p[A, B]), in general (AAAB) ^ (ABAA) unless A and B commute 
or one of them commutes with p. 

7 To show that Tr(pS) = 0, we let B = P^np and use the identity BBB = which fol- 
lows from B 2 = B, BB + BB = B, BBB + BB = BB. Let P a be the one-dimensional 
projectors |a)(c«| onto the eigenvectors of p with non-zero eigenvalues p a (repeated if de- 
generate). Then, P a Pp = S af3 P a , B = J^pPp, P a B = BP a = P a , p = Y. a P aPa> 
Trp = p a = 1, Y2 a p a = 0, S = — k B Pp m P^i always well defined because the 
sum is restricted to the nonzero P/3's, S = — k B PfiVfilVli ~ k B ^2p Pplnpp, and fi- 
nally Tr(p5) = -k B ~ k B E Q E/3 ^{P a Pi3)p a \np p = -k B J2 a Tr(P a P a )p a \np a = 
because TT(P a Pp) = <5 Q(3 Tr(P Q P Q ) and = Tr(BBBP a ) = Tr(BBP a ) = Tr(P a B) = 
E^Tr(P Q P )3 ) = Tr(P a P a ). 
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So, thinking geometrically, 7(f) is a parameterized path in the (square root) 
density matrix space and j(t) gives the tangent vector to this path. If the 
hermitian operator A is not varying directly with time, A'(t) ■'jit) gives the rate 
of change of its mean value (A)(t) as one follows the given path. To maximize 
the rate of change of the mean value, the tangent vector 7(f) should be parallel 
to A'(t), whereas to hold the mean value constant, it should be perpendicular. 
From this follows the interpretation of 2-f, H' and S' as the normalization, 
energy, and entropy gradients, respectively. For a trajectory j(t) to conserve 
normalization and energy, the tangent vector j(t) must be always perpendicular 
to 7 and H'. This holds automatically for unitary evolutions, for which j(t) 
is perpendicular also to S' (see Section |2~3)) . For a more general evolution, to 
maximize the generation of entropy, 7(t) should be parallel to 5", but in general 
this is in conflict with the conservation of normalization and energy. So, in 
Section l2~4l we take j(t) to be parallel to the component of S' perpendicular to 
both 7 and H', i.e., in the direction of greatest entropy increase consistent with 
the constraints of normalization and energy conservation. 

2.3 Hamiltonian Evolution and Time-Energy Uncertainty 
Relations 

Eqs. (fTT|) and (fT5|) imply that 7 remains unit norm and d(H)/dt = (dH/dt) 
(energy conservation when 7 is orthogonal to both 7 and the energy gradient 
vector H'. 

This is the case for purely Hamiltonian evolution, whereby 

j H =i 7 AH/h p H = -i[H,p]/h. (23) 

Note that jn is everywhere orthogonal also to S' and hence also the entropy is 
time invariant. It is actually so in a very special way, as each eigenvalue of p 
is time invariant under unitary evolution. Note also that, since jAH • jAH = 
(AH AH) — h 2 j H ■ j H , the Fischer-Rao metric dt — 2y/-jH ■ -jn dt = dt/rH 
defines an intrinsic Hamiltonian time th such that 

(AH AH) t h = h 2 /A . (24) 

This can be interpreted as the shortest characteristic time of unitary evolution, 
because indeeed the standard (Mandelstam-Tamm-Messiah [45]) time-energy 
uncertainty relation can be stated as follows [44J [46] , 

t 2 fh = (AFAF)/\d(F)/dt\ 2 >ijj (AH AH) r 2 FH > h 2 /A , (25) 

where F is any hermitian operator and tfh is the characteristic time of change 
of its mean value (F) under Hamiltonian evolution. A well-known implication 
of l|25p is that the mean value (F) of any observable cannot change at a rate 
exceeding 2^/ (AFAF) (AH AH) /%. 

8 In particular, (H) is time invariant if H is time independent (isolated system). 
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2.4 Time Evolution Along the Path of Steepest Entropy 
Ascent 

Instead of pure Hamiltonian evolution, let us assume 

7 = in + To , (26) 

with jd in the direction of steepest entropy ascent compatible with the con- 
straints jz> ■ 7 = (conservation of Trp = 1) and jjj ■ H' = (energy conserva- 
tion). Thus, we assume 7r> orthogonal to 7 and if'. As a result it also turn out 
to be orthogonal to jn ■ 

To identify the direction of steepest entropy ascent, we follow a simple ge- 
ometrical construction based on the well-known but seldom used standard ge- 
ometrical notions reviewed in the Appendix, which from here on we give for 
granted. Let L(y,H') denote the real linear span of vectors 7 and H' . Denote 
by S' L the orthogonal projection of the entropy gradient vector S' onto L, and 
by S'^ L its orthogonal complement, so that 

S' = S' L + S' ±L . (27) 

Assume for simplicity that 7 and H' are linearly independent (the case in which 
they are not is easily covered as done in the Appendix). Then, we may write 



"J ±L(7,_ff') — 



S' 


7 


H' 








S' 


■ 7 


7 ' 7 


H' 


7 




(AS)' 


(AH)' 


S' 


H' 


7 -if' 


H' ■ 


H' 




(AS)' ■ (AH)' 


(AH)' ■ (AH)' 




7 


• 7 7 


H' 






(AH)' 


(AH)' 




7 ' 


H' H' 


■ H' 











(28) 



(AS)'.(AHy (AH y_ (A ay (ASAH) 
~ {AS) - (AH)'-(AH)' {AH) - {AS) (AHAH) [AH) ' 



(29) 



Clearly, operator S' ±L is the component of the entropy gradient operator S' 
orthogonal to both 7 and H'. Using the above expression, it can be read- 
ily verified that S' ±L is orthogonal also to jh (|23| . Moreover, S' ±L ■ S' ±L — 
(AS AS) - (AH AS) 2 /(AH AH). ' 
Therefore, we assume that jo is in the "direction" of S' ±L , and we let 

• 1 S'j_ LhMI) 

where td is some positive constant or functional of p which we call the "intrinsic 
characteristic time of dissipation" also for the reasons to be further discussed 
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below. As done in our original work, for simplicity of notation we also define 
the positive characteristic time functional 



T = T Dy S'±L{~f,H') ' 5" ±L(~/,H')/k B , where k B is Boltzmann's constant. (31) 

2.5 Dynamical Equation for the Density Operator of an 
Isolated System 

When "translated" in density operator formalism, our equation of motion (|26p 
with jjj given by (|30p may be written in the following compact form 

p = - l -[H,p] + 1 (AM,/)} , where M = and 8 H = 

h l ' FJ 2T Dy /(AMAM) e H (AS AH) 

(32) 

Operator M, that we call the "non-equilibrium Massieu operator" , is a nonlin- 
ear function of p not only through the logarithmic dependence in operator S but 
also through the nonlinear functional 9h, which we may call "constant-energy 
noncquilibrium temperature" , because it will smoothly approach the equilib- 
rium temperature as the state will approach at constant energy the maximal 
entropy stable equilibrium state (see below). Notice also that any reference to 
the square root operator 7 disappears, proving that its use is only auxiliary to 
the given geometrical construction, as it is auxiliary to the equivalent variational 
formulation given in [37] (see below). For this reason, the first formulation [B] 
avoided the explicit use of ^fp. 

3 Features of the proposed dynamics for a single 
isolated particle 



Detailed proofs of the general features of Eq. (|32|) that we outline in this section 
are given in Refs. EH EQ SH [46] . 



3.1 Conservation of Nonnegativity of the Density Opera- 
tor 

Eq. (|32p generates a smooth continuous path in state space and maintains at 
zero any initially zero eigenvalue of p. As a result, no eigenvalue can cross 
zero and become negative, neither in the future nor in the past, thus assuring 
conservation of the nonnegativity of p both in forward and backward time. 



3.2 Entropy Generation Rate 

The rate of entropy change (more precisely, entropy "generation" , since so far 
Eq. (|32p is meant to be a model for an isolated or an adiabatic system) takes 
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the following various explicit expressions 

d(S) d{-k B Trpbxp) , 1 . 

= ^ = 7 -S = I ^S ±L -S ±L =Ak B r 7D - lD (33) 

^ ^(AMAM) = -L ((A5A5) - <^>) = ^^(AMAM) , 

(34) 

and is clearly positive semi-definite owing to the well-known positive semi- 
definitcness of Gram determinants (Appendix) and scalar product norms [see 
the last two of Eqs. 



3.3 Characteristic Times and Time-Entropy Uncertainty 
Relation 

Because of the orthogonality between 7£> and jn, if at one time [H, p] = then 
p(t) commutes with H at all times. For such particular, "purely dissipative" 
trajectories, the Fisher-Rao metric takes on the following interesting explicit 
expressions 



= 2,-T^&=^m a= ^m dt (35) 




yj (AM AM) A _ 1 / /A0AC , X (AH AH) ^ _ dt 

_ 



k B r dt -l-r^ S ^--^ dt --^ ^ 



where the last equality justifies our calling tjj the natural "intrinsic dissipative 
time" . The last of Eqs. ([35]) shows that the alternate characteristic time r is 
directly related to the rate dl/dt at which the state operator 7 follows the steep- 
est entropy ascent trajectory, modulated by the dimensionless entropy gradient 
d(S/k B )/d£ along the trajectory, 

1 di/dt 

T d(S/k B )/dl ' 1 ' 

so, we see that when time intervals are measured in units of the "local" (or 
instantaneous) value of r and entropy in units of k B , the "speed" along the 
steepest entropy ascent trajectory (geodesic) coincides with the local slope of 
the entropy surface along the trajectory, 

de _d(s/k B ) 



dtjr dl 

which again justifies the interpretation of r as an intrinsic dynamical time. 
Finally, we notice that dSj dl equals half of the norm of the component of the 
entropy gradient operator 5" orthogonal to the linear manifold L(j, H') defined 
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by the gradients of the constraints (here, for simplicity, only normalization and 
energy conservation), 



dS/de = J yJS' ±L ^, HI) ■ S' ±L(1>W) . (39) 

A noteworthy result follows from Eq. (|34f) together with the general inequal- 
ity (AS AS) > (AM AM) (Appendix) . By defining the characteristic time ts of 
the entropy generation, we find the following general time-entropy uncertainty 
relations 

, (AS AS) (AM AM) {k B rf 2 
Ts ~ \d(S)/dt\i ~ \d{S)JW = AM) ~ (AM AM) t s > (k B r) 

(40) 

and also 

(AS AS) r| > (k B T) 2 , (41) 

which imply that the rate of entropy generation cannot exceed the following 
bounds, 

d(S)/dt < ^/(ASAS)/t d < (ASAS)/k B r . (42) 



3.4 Equilibrium States and Limit Cycles 

From Eq. p3|) we see that the rate of entropy generation is zero (and the 
evolution is Schrodinger-von Neumann) if and only if 70 = 0, i.e., when 5" lies 
in L(7, H'). Then, the density operator may be written as 

B exp(— g/ k B T)B 2 
P= Tr[BeM-H/k B T)B] > iOTSOmCB = B ' (43) 

where T — H = (AH AH) / (AH AS) =6 S = (AH AS) / (AS AS) = ^ (AH AH) /(AS AS). 
We call these the "nondissipative states" . Proofs of the above and following re- 
sults are straightforward, and detailed in the original papers [9l HU1 HT1 130] . 

Because 70 = 0, nondissipative states remain nondissipative at all times, 
— 00 < t < 00. Therefore they obey unitary Hamiltonian evolution. If [B,H] = 
they are equilibrium states. If [B,H] ^ their unitary evolutions are limit 
cycles of the dynamics, and B(t) = U(t)B{0)U- 1 (t) with U(t) = exp(-itH/h). 
Limit cycles can be mixed, if Tr£> > 1, or pure, if TiB = 1. The latter case 
coincides with the usual Schrodinger dynamics of standard quantum mechanics. 
Except for when B = I (the identity operator) all these equilibrium states and 
limit cycles are unstable (according to Lyapunov). 



3.5 Dynamical Group, Not a Semi-Group 

As proved in the original papers and again in Refs. [T71GI7], the solutions of our 
nonlinear dynamical equation form a group, not a semi-group, of "trajectories" 
in density operator space. This is so because of a very important feature of 
the dynamics, namely, that along every trajectory the zero eigenvalues of p are 
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invariant, the range of p and its cardinality TrPR anp are invariant. However, 
the nonzero eigenvalues of p get smoothly rearranged towards higher entropy 
distributions approaching partially canonical distributions of the form 1431 Be- 
cause positive eigenvalues of p remain positive at all times, it follows that every 
trajectory is unique and well-defined not only in forward time but also back- 
wards in time (an explicit proof of the uniqueness feature is given in [17] for a 
two level system, and is also discussed in general in [37]). 

We therefore have a "strongly causal" group of dynamical evolutions, with 
inverse defined everywhere, unique trajectories through every state p, fully de- 
fined both forward and backward in time, thus allowing full reconstruction of 
the past from the present. It is an example of an invertible dynamics which nev- 
ertheless is fully compatible with all thermodynamics principles and in addition 
is largely irreversible in that in forward time and for a system which is isolated 
(or "adiabatic" in the sense that the Hamiltonian operator may be time depen- 
dent) entails and describes explicitly entropy generation along the direction of 
maximal increase. 

A remark about invertibility is in order. An often encountered, misleading 
assertion is that to be "irreversible" a dynamical evolution must not have an in- 
verse. To justify the theory of dynamical semi-groups, the misleading suggestion 
has often been made that non-invertibility is an indispensable feature for the 
description of thermodynamic irreversibility, so that the equation of motion can 
be solved only forward in time, not backwards, and causality is thus retained 
only in a "weak form": future states can be predicted from the present state, 
but the past cannot be reconstructed from the present. 

Our nonlinear dynamical group challenges this idea. The existence of ther- 
modynamic irreversibility is not incompatible with causality in the strong sense: 
it does not necessarily forbid the possibility to reconstruct the past from the 
present. If p(0) and p(t) are the states at times and t, they are related by the 
dynamical map p(t) = A t (p(0)) through the solution of the equation of motion 
for the time interval from to t with initial condition p(0). If the inverse map 
exists, it points from the final state back to the initial state, p(0) — A^ 1 (p(t)), 
and there is a one-to-one correspondence between initial and final state. This 
is our case. It shows that existence of the inverse map, does not rule out the 
possibility that the functional which represents thermodynamic entropy could 
be non- decreasing in forward time. In Ref. [17] we present some numerical 
solutions which exemplify how any given distribution of eigenvalues belongs to 
a unique smooth solution which among other features identifies (as t — > — oo) a 
lowest-entropy (not necessarily zero-entropy) "ancestral" or "primordial" state. 

3.6 Stability of Equilibrium States and the Second Law of 
Thermodynamics 

Each partially canonical equilibrium density operator of the form (|43[) maxi- 
mizes the entropy when restricted to the subset of density operators that share 
the same kernel. Such states are equilibrium (i.e., time invariant if H is time 
invariant), but are unstable whenever at least one eigenvalue of p is equal to 
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zero. In fact, a minor perturbation which changes the zero eigenvalue to an 
arbitrarily small nonzero value, would proceed away towards a quite different 
equilibrium of higher entropy. Instead, any trajectory with no null eigenval- 
ues, maintains such feature at all times (invariance of the cardinality of the 
set of eigenvalues), and approaches in forward time the unique, fully canonical, 
maximum-entropy density operator compatible with the initial values of the 
constraints (remember that for simplicity we are considering here only systems 
whose approach to equilibrium is constrained only by energy conservation and, 
of course, normalization). 

Therefore, the only dynamically stable equilibrium states (again, stable ac- 
cording to Lyapunov) are those given by Eq. (|43[) with B = I. There is only one 
such canonical density operator for every value of the mean energy (H), which 
through Tr Hp — (H) fixes the temperature T in (|43[) . Existence and uniqueness 
of stable equilibrium states for every value of the energy is the essence of the 
Hatsopoulos-Keenan statement of the second law of thermodynamics, which we 
may state as follows [151 1401 148) : Among all the states of a system that have a 
given value of the energy and are compatible with a given set of values of the 
amounts of constituents and the parameters of the Hamiltonian, there exists one 
and only one stable equilibrium state. From this statement of the second law, the 
Kelvin-Planck, the Clausius, and the Caratheodory statements can all be shown 
to follow as logical consequences (explicit proofs of this assertion can be found in 
[IH p.64-65 (Kelvin-Planck), p.133-136 (Clausius), and p.121 (Caratheodory)]). 

This statement of the second law brings out very clearly the apparent conflict 
between mechanics and thermodynamics, a contrast that for over a century has 
been perceived as paradoxical. In fact, within mechanics, classical or quantum, 
the following so-called minimum energy principle applies: Among all the states 
of a system that are compatible with a given set of values of the amounts of 
constituents and the parameters of the Hamiltonian, there exists one and only 
one stable equilibrium state, that of minimal energy. Comparing this assertion 
with the statement of the second law just reviewed, leads to a paradox if we in- 
sist that the two theories of Nature contemplate the same set of states. Indeed, 
for fixed amounts of constituents and parameters of the Hamiltonian, mechan- 
ics asserts the existence of a unique stable equilibrium state (that of minimal 
energy), whereas thermodynamics asserts the existence of infinite stable equi- 
librium states (one for every value the mean energy can take). 

The paradox is removed if we admit that the "pure" states contemplated 
by Quantum Mechanics are only a subset of those contemplated by Thermody- 
namics. This resolving assumption was very controversial when Hatsopoulos and 
Gyftopoulos first introduced it in [15] . However today — more as a byproduct of 
the more recent vast literature on quantum entanglement and quantum informa- 
tion than as a result of thermodynamic reasoning — an assumption to this effect 
is ever more often being included in the postulates of quantum theory (compare 
for example the postulates of quantum theory as stated, e.g., in the recent m] 
with those stated in 1968 by Park and Margenau [50]). The discussions on the 
relations between this fundamental assumption and thermodynamics is flourish- 
ing in the physics literature (see [HI [51] and references therein) . Unfortunately, 
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pioneering contributions such as [TS] are seldom acknowledged. 

As mentioned in the introduction, the validity of the present steepest-entropy- 
ascent or maximal-entropy-generation mathematical formalism even outside of 
the original framework for which it was developed [131 EE M US] has been recog- 
nized and suggested by this author long ago, not only for quantum dynamical 
phenomenological modeling [3U] but also as a general tool for modeling relax- 
ation and redistribution of nonequilibrium probability or other positive-valued 
distributions in a variety of fields [IHl HH] • 

3.7 Variational Formulation 

In 2001, Gheorghiu-Svirschevski [37] re-derived Eq. ([52"]) from a variational 
principle that in our notation is 

d(S) , . d(H) n dTrp n 2 

max — - — subject to — - — = 0, — - — = 0, and ■yo • Id = c , (44) 
dt dt dt 

where the last constraint signifies that we maximize the rate of entropy gen- 
eration at fixed norm of the operator 70, hence we are free to vary only its 
direction (c 2 is some real functional independent of 70). Introducing Lagrange 
multipliers, 

L = j D ■ S' - Ai j D ■ 7 - X H j D ■ H ' - A r 7 D • j D . (45) 

and maximizing L with respect to 70 yields 

S' - Ai 7 - \h H' - 2A T 7d = , (46) 

where the multipliers must be determined by substitution in the constraint 
equations. It is easy to verify that our expression of 70 in Eqs. ([28]) and (|30|) 
yields the explicit solution of Eq. ([46)) . Using ([36)) we see that with c 2 = 1/4t|) 
we get exactly our quantum dynamical evolution equation. 



3.8 Onsager Reciprocal Relations even Far from Equilib- 
rium 

Any nonequilibrium p can be written as 

p = iK t±lll 3 ' (47) 

9 TrBeM-EjfjXj) ' 1 ' 

where the set {I,Xj} spans the real space of hermitian operators on H, and 
B = B 2 is a projector (actually, B = Pjun P )^ We may call the set {X,} a 
"quorum" of observables, because the measurement of their mean values {(Xj)} 



9 To prove Eq. I|47| l. let Pk B tp denote the projector onto the kernel of p, i.e., the eigenspace 
belonging to the zero eigenvalue (if p is non-singular, PKerP projects onto the null vector of 
H), and let -B('H) denote the real space of the hermitian operators on 7i, equipped with the 
inner product X -Y = Tr(XY) and let {/, Xj} be a basis for B(7i). The (nonnegative definite) 
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fully determines the density operator. The empirical determination of a quan- 
tum state has been recently called "quantum tomography" . In an almost for- 
gotten seminal series of papers in 1970-1971, Park and Band [52] devised elegant 
systematic rules to construct such a quorum of observables. 
Given such a quorum, we can write 



(Xj) = TripXj), (S) = k B f + k B ^ f 3 (X,), and k B f, 



d(Xj) 



, (48) 



where k B fj may be interpreted as the "generalized affinity" or force, conjugated 
with the observable associated with operator Xj. 

Let us focus on the dissipative term in our equation of motion (|32[) and the 
rate of change it induces on the mean value (Xj) of each quorum observable. 
We call it the "dissipative (part of the) rate of change" of observable Xj , 

(Xj) D = - lD ■ X' 3 , with X' 3 = 2 1 X J . (49) 

From the expressions we derived for 7u in our steepest entropy ascent dynamics, 
we find the following linear relations between dissipative rates and affinities 

(X i ) D = J2fjL ij (p) , (50) 

3 

where the coefficients are nonlinear functionals of p which form a symmetric, 
positive semi-definite Gram matrix [{L^ (/?)}], "generalized conductivity", 



Lij(p) 



{AXiAXj} (AHAXj) 
(AX.AH) (AHAH) 



t (AHAH) 



- Lji(p) , (51) 



As a result, the entropy generation rate may be written as a quadratic form in 
the affinities 

i 3 

When [{Lij(p)}} is positive definite, we denote its inverse, "generalized resis- 
tance", by [{Rij(p)}], 

/i = EWi)i): (53) 

i 

and the entropy generation rate may then be written also as a quadratic form 
in the dissipative rates 



operator — ln(p + PKerp) is well defined for every p and, since it belongs to P(W), we may 
write it as - ln(p + P Kc r P ) = hi + J2j fj x j- Therefore, p + P Kc r P = e~ fo exp(- J2j fj x j)- 
Multiplying by P Ranp and using the identities PR anp P K crp = and P Ranp p = pPR anp = P, 
we obtain p = e~ ^°PR anp exp(— ■ fjXj). Finally, by imposing Trp = 1, we find e'° = 
Tr[P Ranp cxp(-E,/i^')]- 



20 



as well as as a sum of the dissipative rates of change of the quorum observables 
each multiplied by its conjugated affinity 

^ = k B J2fi{Xi)n. (55) 

i 

Notice that the parametrization of density operators given by Eq. [47] in 
terms of the real variables fj, implies that we can write the entropy as (S) = 
(S)(fi, . . . , /n 2 -i) m view of the fact that our dynamics conserves the cardinal- 
ity of p (and hence TtB is invariant). 

3.9 Nonlinear Master Equation for Energy Level Occupa- 
tion Probabilities 

When written for the nm-th matrix element of p with respect to an eigenbasis 
{\e 3 }} of H, Eq. ((32j) becomes 

dgnm i !T? „ s . 1 * ( a Ae„ + Ae m N 

= -^Pnm[E n -E m ) + y^2-^ u nrU mr p r ^As r — j , (56) 

where Ujk = (^j\Vk), {\Vk}} is an eigenbasis of p, p^s its eigenvalues, Ask = 
Sk - (S), Sfc = -k B liip k if p k ^ 0, s k — if p k = 0, ej the eigenvalues of H, 
Aej = ej — (H), and Oh , k B and r as already defined above. From (|56|) we see 
that if at one instant of time, [H, p] = and we select a common eigenbasis, 
then Ujk = Sjk and dp nm /dt = for w/ra, which means that the condition 
[H, p] = holds along the entire trajectory. In such special but nontrivial cases, 
the eigenvalues of p get redistributed according to the nonlinear master equation 

dp n 1 ( * Ae„\ 

-J!" = TT^P" As ™ - -5— . ( 57 ) 



dt fc B r \ 6h 

whose fundamental features are analyzed and numerically exemplified in Ref. 
|17j , where we wrote it in the following equivalent form 

"TT = - ~ \Pn m Pn + ap n + P e n p n ] , (58) 
at t 

with the nonlinear functionals a and (3 defined by 

_ E» apt Ej- e jPj ln Pj ^ E» gj m P» Ej e jP 3 _ EiP» Ej e jP 3 - E» e »p» ln P» 
Ei e iP« - ( Ei e ^J E 4 e fpi - ( E< e iP 

(59) 

In the usual statistical mechanics framework, the eigenvalues of p when 
[p, H] = are interpreted as "occupation probabilities" , meaning that the sys- 
tem is to be thought of as in a particular, unknown, pure state; then, these 
"probabilities" are understood as an expression of the uncertainty as of which 
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pure state the system is actually "occupying". In our original framework, in- 
stead, the eigenvalues of p when [p, H] = are interpreted as "degrees of energy 
load sharing" among the different modes (eigenvectors of H) with which the 
system can internally accommodate its mean energy. From this point of view, 
entropy measures an overall degree of sharing between the available and ac- 
tive modes (i.e., those with non-zero eigenvalues). Entropy generation measures 
therefore the rate at which the spontaneous internal dynamics redistributes en- 
ergy among the available modes, to achieve maximal sharing. 



4 Extension of Equation (1321) to Composite Sys- 
tems 

The nonlinear, dissipative term jjj in Eq. (126|) provides a strong coupling be- 
tween the energy storage modes of the single-particle system, additional to the 
coupling entailed by the linear, unitary term 7^ through the structure of the 
particle's Hamiltonian operator H. Were we to apply Eq. (f2T)]) without modifi- 
cations to a system composed of two particles A and B (or to a more complex 
composite system) the term would couple the subsystems and make them 
exchange energy even in the absence of an interaction term in the Hamiltonian 
H , thus violating both separability and no-signaling criteria. Because of the 
nonlincarity which is intrinsic in the steepest entropy ascent construction, if the 
model equation of motion is to meet these criteria, for a composite system, the 
structure of the interactions and the internal constraints between subsystems 
must be described not only through the Hamiltonian operator, but also through 
the structure of the dynamical equation itselfF°l 

Suppose Alice and Bob, A and B, are the two elementary subsystems of an 
adiabatic system. Each subsystem is a single particle. Alice and Bob may be 



interacting 
nonintcracting 
correlated / entangled 
uncorrelated 



H j£ H A ® Ib + Ia<E)H b ; 
H = H A ®I B +I A ®H B \ 
S(p) ^ S(p A ) ®I b + Ia® S(p B ); 
S(p) = S(p A ) ®I b + Ia® S(p B ); 



10 It is noteworthy that if the proposed nonlinear evolution law is supposed to be a fun- 
damental law of nature (i.e., not just a phenomenological modeling tool), then one should 
specify criteria for dividing a system into its separate elementary constituents. A unitary 
Hamiltonian dynamics depends only on the Hamiltonian operator H, regardless of the level 
of description, i.e., of whether we reach the given H by considering as elementary con- 
stituents the individual atoms, or the individual electrons within an atom and the nucleus, 
or the nucleons, or quark and gluons, or the electron field and quark field, etc. For the 
given H and a given mean value (H) of the energy, also the unique stable equilibrium state, 
P{H) = ex P( — P(H) H) /Tr exp(— P(B)H), is independent of the level of description. But if the 
relaxation to stable equilibrium is described by our nonlinear law, then the dynamics depends 
strongly on the assumed level of description, because through a unique internal relaxation 
time for each elementary subsystem, this dynamics fully couples all the internal modes of 
the elementary subsystem in a local effort to follow a path of steepest ascent in the locally 
perceived value of the overall entropy (see below). 
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where here S(p) denotes as before the operator — k B PR anp Inp, p is the density 
operator of the overall system, and pa, Pb the reduced local density operators. 

Our construction [TO] was designed so as to obtain a dynamical system 
obeying the following separability and no-signaling criteria [16] . 

• For permanently non-interacting subsystems A and B, every trajectory 
passing through a state in which the subsystems are in independent states 
{p = Pa ® Pb, where pA = Tr^p and ps = Tr^p) must proceed through 
independent states along the entire trajectory, i.e., when two uncorrelated 
systems do not interact with each other, each must evolve in time inde- 
pendently of the other. 

• If at some instant of time two subsystems A and B, not necessarily non- 
interacting, are in independent states, then the instantaneous rates of 
change of the subsystem's entropies — k B Ti(pA lnp^) and — k B Ti(pB hips) 
must both be nondecreasing in time. 

• Two non- interacting subsystems A and B initially in correlated and/or 
entangled states (possibly due to a previous interaction that has then 
been turned off) should in general proceed in time towards less correlated 
and entangled states (impossibility of spontaneous creation of any kind of 
correlations) . 

• When subsystems A and B are not interacting, even if they are in entan- 
gled or correlated states, it must be impossible that the time dependence 
of any local observable of one subsystem be influenced by any feature of 
the time evolution of the other subsystem (no-signaling condition) . 

Notice that we do not request that existing entanglement and/or correla- 
tions between A and B established by past interactions should have no influence 
whatsoever on the time evolution of the local observables of either A 01 B. In 
particular, there is no physical reason to request (as is often done) that two 
different states p and p' such that p' A = pa should evolve with identical local 
dynamics (dp' A /dt = dpA/dt) whenever A does not interact with B, even if en- 
tanglement and/or correlations in state p differ from those in state p' . Rather, 
we see no reasons why the two local evolutions could not be different until 
spontaneous decoherence (if any) will have fully erased memory of the entan- 
glement and the correlations established by the past interactions now turned 
off. In fact, this may be a possible experimental scheme to detect spontaneous 
decoherence. In other words, we will not assume that the local evolutions be 
necessarily Markovian. 

Compatibility with the predictions of quantum mechanics about the gener- 
ation of entanglement between interacting subsystems that emerge through the 
Schrodinger- von Neumann term —i[H,p]/h, requires that the dissipative term 
may entail spontaneous loss of entanglement and loss of correlations between 
subsystems, but should not be able to create them. 
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To this end, we devised [5] OH] a construction which hinges on the definitions 
of the following "locally perceived energy" and "locally perceived entropy" op- 
erators, as a result of which our composite dynamics implements the ansatz of 
"steepest locally perceived entropy ascent" , 

(AH) A = Tr B [(I A <g> p B )AH] and {AH) B = Tr A [(p A <g> I B )AH] (60) 
{AS) A = Tr B [(I A ®p B )AS] and (AS) B = Tr A [[p A «8> I B )AS] . (61) 

A geometrical construction (details in [5J [TOl HI] ) analogous to that outlined 
in Section[2]for a single particle, leads us to a composite-system steepest-locally- 
perceived-entropy-ascent dynamics with the form 

^ = -l[H, p] + —r- — {(AM) A , PA }® PB + ^—PA ® {{AM) B , PB } , (62) 
dt n 2k b t a 2k b t b 

where t Ai t b are local characteristic times and, for J = A, B, 

(AM) J = (AS)' 1 - (AH) J /6 H , with HJ = {(AH) J (AH) J ) / {{AS) J (AH) J ) . 

(63) 

Each local dissipative term separately "conserves" the overall system's mean 
energy (H) = Tr(pH). Each subsystem's contribution to the overall system's 
rate of entropy change is positive semidcfinitc 

^ = —J— ((AM) A (AM) A ) + -1— ((AM) B (AM) B ) . (64) 
at k B r A k B T B 

If Alice and Bob interact, it is only the Hamiltonian term in the evolution 
equation which during the interaction builds up correlations. Once generated, 
these correlations survive even after A and B separate, even if the loose touch 
completely When that happens, A and B remain correlated but begin to evolve 
independently of one another. This is reflected in the local structure of our 
equation and in particular of operators (AS) A , (AS) B , (AH) A , (AH) B . 

Despite the nonlincarity, the equation prevents no-signaling violations, in 
that it satisfies the following strong separability conditions. Namely, denoting 
by p AB (p, H) the rhs of the equation, it is easy to show that, for any p and any 
H A , H B , 

Tr B [p AB {p,H A ®I B +I A ®H B )} = f A {{AS) A ,H A ) , (65) 
Ti A [ PAB (p, H A ®I B + I A ®H B )} = f B ((AS) B ,H B ) , (66) 



Conditions (|65j) and (|66|) . when restricted to uncorrected states, p — p A ®p Bl 
define the conditions of weak separability ; which of course are a corollary of 
strong separability. 

However, existing correlations do influence the local evolutions, which there- 
fore are not Markovian in that they do not depend only on the respective local 
(reduced) states p A and p B . 
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5 Additional Phenomenological Modeling Equa- 
tions for the Density Operator of a Coupled 
System 

In this section, we return to the problem of describing the effective interaction 
between a system and a reservoir. But instead of starting from Hamiltonian 
dynamics and adopting suitable approximations so as to arrive at the KSGL 
equation as discussed in Section ITTTT we take a fully phenomenological approach. 
Using our geometrical construction, it is easy to "design" dynamical equations 
that exhibit dynamical features that we expect from typically thermodynamical 
energy balance and entropy balance considerations. So, we may obtain a variety 
of similar dynamical equations that, though perhaps not as fundamental, may 
nevertheless be very useful in the phenomenological description of nonequilib- 
rium phenomena of non-isolated systems. 

For example, recent major advances in micro- and nano-technological ap- 
plications, often call for a detailed description of the time evolution of non- 
equilibrium states that are far from thermodynamics equilibrium and cannot be 
described by partially canonical entropy density operators. In such far noncqui- 
librium regime, the assumption of linearity underlying the standard theory of 
irreversible processes may be cease to hold. Yet we may need to describe the 
simultaneous energy and entropy exchange which occurs between our quantum 
system and, say, a thermal reservoir at temperature Tq, whereby the ratio of 
the energy to the entropy exchanged is equal to Tq. 

Thus, we will consider additional terms of the form 

1 {AG, p} , with G = S-^r , (67) 



2t g ^(AGAG) 

where the dots represent other terms as in (p?2"j) or even in the KGSL Eq. (fl]), 
G is another " non-equilibrium Massieu operator" that depends on the choice 
of functional 9 (see below), and tq is the characteristic time of decrease of (G). 
For 6 = 6h, G coincides with operator M in Eq. (p{2")) . Notice that if 9 is chosen 
to be a constant, then T G d(G) /dt = —\J {AG AG) < except for jG' = 0, that 
is, jS' = jH'/O or (AH AH) = 9 2 (ASAS) which occurs only for states of the 
form g3]) with T = 9. 

Eq. (|67p generates energy and entropy rates of change according to 



d(H) (AH AH) ( 1 \\ , d(S) (AS AS) , 

- 1 1 and — ; — = - 1 — 



dt r Gv /(AGAG) \9 H 0J dt t g ^/(AGAG) 

(68) 

where 9 H = (AH AH) /(AH AS) and 9 S = (AH AS) /(AS AS) as defined be- 
fore. Note that in general 9^ > 9g with strict equality only at states of form 
P5|) . and 9tf9g > which means they always have the same sign (that of 
(AH AS)). 
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Alternatively, we may also consider equations of the form 

1 



2r FV /(AFAF) 



{ AF, p} , with F = H -OS , (69) 



where F is a "non-equilibrium Helmholtz free energy operator" that depends 
on the choice of functional 0. In this case, 

d(H) (AH AH) / 6 \ , d(S) (9 s - 6) (AS AS) 

1 — - — and 



dt T p \J (AFAF) V <W * Tp^J (AF AF) 

(70) 

Eqs. ({69]) and {70j) are related to (67]) and (J6SJ) by the fact that (AFAF)f(AGAG) 
(AHAH)f(ASAS) = 6> ff (9 s . 

5.1 Smooth Isoentropic Extraction of the Adiabatic Avail- 
ability 

When a system is in a noncquilibrium state, we call "adiabatic availability" the 
largest amount of energy that can be extracted in the form of work without leav- 
ing any other effects external to the system and without changing the system's 
Hamiltonian operator^ It is given by (*) = (H) - (H) s where (H) = Tr[Hp] 
the mean energy of the nonequilibrium state and (H) s = Tr[H p s (H)] the mean 
energy of the unique stable equilibrium state p s (H) — exp(—(3 s H) /Tr exp(— (3 S H ) 
that has the same entropy as the given state p, i.e., such that (S) — (S) s . As 
thoroughly discussed in [40] in general terms, and in [15l [53] in the quantum 
framework, the adiabatic availability cannot in general be completely extracted 
by means of a unitary evolution, owing to the fact that a unitary process cannot 
change the eigenvalues of p. Instead, a process is required that while maintaining 
the entropy invariant, smoothly modifies the eigenvalues of p until they become 
canonically distributed. At the end of this isoentropic change of state, p has the 
form (43]) with B = 1. 

Whereas finding a practical way to control and interact with the system's 
dynamics so as to extract its adiabatic availability from an arbitrary initial state 
may be a very hard problem, our geometrical construction makes it straight- 
forward to design a dynamical equation that describes phenomenologically such 
an extraction, along a steepest-energy-descent trajectory at constant entropy. 
It suffices to take 7 proportional to —H' ±L ,^ s ,s, i.e., the component the energy 
gradient H' orthogonal to both 7 and S', 

{AH)' (AS)' 
(AH)' ■ (AS)' (AS)' ■ (AS)' 



H'±L(y,S') 



-(AHY {AHAS) (A^Y 
(AS)' ■ (AS)' 1 ' (AS AS) K ' ' 

(71) 



11 Even if the Hamiltonian operator is a function H(X) of some controllable parameters, so 
that the state of the system is given by (rho, A), the adiabatic availability is the largest work 
that can be extracted with no net changes in A. 
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In terms of Eq. (|69p for the density operator, we may describe this by choosing 



= 9 S 



(AHAS) 
(ASAS) 



F = = H -9 S S and r F = r K 



(72) 



therefore we may call 9s the "constant-entropy nonequilibrium temperature" 
and F^ the "constant-entropy nonequilibrium Helmholtz free energy operator" . 
With this choice of 9 in Eq. (|69|) . the entropy remains constant while 9s and 
9h smoothly approach the temperature of a canonical or partially canonical 
final state of lowest energy for the given entropy. Notice that the rate of energy 
change d(H)/dt is negative semidefinite (even for states with negative 9h)- 
Because the zero eigenvalues are time invariant here like for Eq. (|32|) . this term 
will extract the full adiabatic availability only if the initial p in non-singular. 

5.2 Smooth Extraction of the Available Energy with Re- 
spect to a Reservoir 

When a system is in a nonequilibrium state or in any state not of mutual equi- 
librium with a given reservoir with temperature Tr, we call "available energy 
with respect to a reservoir with temperature Tr" the largest amount of energy 
that can be extracted in the form of work without any other effects external to 
the combination of the system and the reservoir. It is given by (f2 fl ) = (H) — 
(H)r — Tr ({S) — (S)r) where (H)r and (S)r are the energy and the entropy 
of the unique stable equilibrium state pa = exp(—H/k B Tn)/Trexp(—H/k B Tji) 
with temperature Tr. Again, its definition is discussed in [40) in general terms, 
and in [15] in the quantum framework. 

We can design a dynamical equation that generates a trajectory along a 
smooth descent in available energy by taking Eq. (|69p with 



where Tr is the constant temperature of the reservoir and Fqr is yet another " 
nonequilibrium Helmholtz free energy" . With this choice of 9, the signs of the 
energy and entropy rates (|70|) depend on those of 1 — Tr/0r and 9s — Tr. This 
is a model of a reversible weight process for the system-reservoir composite [40] 
where, by the energy and entropy balance equations, Tad(S)/dt and d(S)/dt 
equal respectively the net rates of energy and entropy exchange (from the reser- 
voir to the system, if positive, from the system to the reservoir, if negative). 
Therefore, the remaining power, —d(H)/dt + Tud(S)/dt = —d(Q R )/dt, is the 
rate of energy extraction in the form of work. From Eqs. ([70]) we see that the 
energy and entropy rates are both zero only at state pr, where 9h = 9s = Tr. 

5.3 Nonequilibrium Heat Interaction 

As a final example, we consider the model of an interaction between our system 
in a nonequilibrium state p and some reservoir (heat bath) at Tq, whereby 
the ratio of the energy and the entropy exchange rates is equal to Tq. The 



9 — Tr ; , F = Fqr 



H — TrS and tr = tr 



(73) 
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usual definition of a heat interaction at Tq (see [3D] for a rigorous definition) 
requires both interacting bodies to be in states very close to their respective 
stable equilibrium states with temperature Tq, because only then the ratio of 
energy to entropy exchanged is equal to Tq. Therefore, the interaction we 
are modeling here is an extension of the standard notion to when one of the 
interacting systems is far from thermodynamic equilibrium (where temperature 
is not defined). It is easy to verify that by taking Eq. (|69|) with 

e = 0Q = 8s "~l Q , F = F Tq =H- 9qS and t f = t f , (74) 

t) S - 1q « 

we obtain a smooth trajectory where at all times d(H)/dt — Tq d(S)/dt. 



6 Conclusions 

In this paper we discuss the geometrical construction and the main mathemat- 
ical features of the maximum-entropy-production/steepest-entropy-ascent non- 
linear evolution equation proposed long ago by this author in the framework of a 
fully quantum theory of irreversibility and thermodynamics for a single isolated 
or adiabatic particle, qubit, or qudit. The same mathematics has been recently 
rediscovered by other authors, with various physical interpretations. 

The nonlinear equation generates a dynamical group, not just a semigroup 
like for KSGL dynamics. It provides a deterministic description of irreversible 
conservative relaxation towards equilibrium from an arbitrary initial density 
operator. It satisfies a very restrictive stability requirement equivalent to the 
Hatsopoulos-Keenan statement of the second law of thermodynamics which 
therefore emerges as a general theorem of the dynamics. It has smooth unique 
solutions both forward and backwards in time. Except for fully characterized 
families of limit cycles and of equilibrium states the entropy functional is strictly 
increasing in forward time and strictly decreasing in backward time. Viewed as 
a model of the relaxation to equilibrium of an isolated single particle system, 
this dynamics entails thermodynamic irreversibility at the single particle level. 

For a multipartite isolated or adiabatic system, we introduce a nonlinear 
projection defining local operators that we interpret as "local perceptions" of 
the overall system's energy and entropy. Each component particle contributes 
an independent local tendency along the direction of steepest increase of the 
locally perceived entropy at constant locally perceived energy. It conserves 
both the locally-perceived energies and the overall energy, and meets strong 
separability and non-signaling conditions, even though the local evolutions are 
not independent of existing correlations. 

In addition, we also show how the geometrical construction can readily lead 
to a variety of thcrmodynamically relevant models, such as the phenomenolog- 
ical descriptions of nonunitary isoentropic evolutions achieving full extraction 
of a system's adiabatic availability or available energy with respect to a reser- 
voir, or the phenomcnological descriptions of a nonunitary nonequilibrium heat 
interaction. 
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Appendix. Orthogonal Decomposition of a Vec- 
tor with respect to a Linear Manifold 

In the paper, we make extensive use of the notation and relations discussed in 
this appendix (see [THUS] for a general but non-quantal context). 
Given a set of vectors go, gi, . . . , gn, the symbol 



L (go, gi, ■ ■ ■ , gn) 



(75) 



will denote their linear span, i.e., the linear manifold containing all the vectors 
that are (real) linear combinations of go, gi, . . . , g n . Given another vector b, 
the symbol 

b L(go,gi,...,g„) ( 76 ) 

will denote the orthogonal projection of b onto the linear manifold L (g , gi, . . . , g„), 
namely, the unique vector in L (g , gi, . . . , g„) such that its dot product with 
any other vector g in L (go, gi, . . . , g„) equals the dot product of b with g, 
i.e., 

g-bi ( g D , gl ,...,g n) = gb (77) 

for every g in L (go, gi, . . . , g„). 

In terms of a set of linearly independent vectors hi, . . . , h r spanning the 
manifold L (go, gi, g n ), where clearly r < n, we can write two equiva- 
lent explicit expressions for the projection {h) L , gQ gi g , of vector b onto 
£(go, gi, ■ ■ ■ , gn)- The first is 



fc— 1 m— 1 



bi(g , gl ,...,g n ) =^X)(b-h fe ) | M ( h i> '-A) 
where M(hi, . . 

M(hi, K) 



(78) 



h r ) is the inverse of the Gram matrix 
hi • hi • • • h r • hi 



hi • h r 



h r • h r 



(79) 
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The second expression is a ratio of two determinants 



J L(go, gi, g„) 






hi 




h r 


b 


hi 


hi - hi 




h r • 


hi 


b 


h r 


hi • h r 




h r • 






hi 


•hi ••■ 


h r 


hi 






hi 


• h r 


h r 


• h r 





(80) 



where the determinant at the denominator, also given by det M (hi, . . . , h r ), is 
always strictly positive because the vectors hi, . . . , h r are linearly independent. 
In the paper, our rate equations are expressed in terms of vectors of the form 



1 


b 


hi 




h,. 


b 


hi 


hi hi 




K- 


hi 


b 


h, 


hi • h r 




h r ■ 






hi 


■hx ••• 


h r 


hi 






hi 


• h r 


h r 


• h r 





b -L£(go,gi,...,g„) - b - b L (g 0igl ,._g n ) 



(81) 

where in writing Equation lSll we make use of Equation lSOl The vector bxL(g , gl gn ) 
is orthogonal to manifold L (g , gi, . . . , g n ); indeed, the vector represented by 
Equation [HI] has the relevant property 

gfc ' b -LL( g0 ,gi,...,g„) = fc = 0, 1, ...,n (82) 

which follows directly from Relation I77i and hence the relation 

b = bi( gD , gll ...,g„) +bxL( g0 , gl ,...,g„) (83) 

represents the unique orthogonal decomposition of vector b with respect to 
manifold L (go, gl, gn)- 

Moreover, we have the other obvious, but relevant properties 

b • bxi(g , gl ,...,g n ) = b xi(g ,g 1 ,..., gn ) • b X z,(g ,gi,...,g„) > > ( 84 ) 

where the strict inequality applies whenever b is not in L (go, gi, . . . , g n ), and 
for any a and b 



a ' b±i(g D)gl) ...,g n ) - a -L£(g 0l gi,...,g«) ' b_LL(g ,gi,...,g„) > 



(85) 



An important formula which derives from Eq. [51] and the usual properties 
of determinants, is 



b±L(g , gi, .... gn) ' bj_L(g , gi, g„) 



det M(b,hi, . . . ,h r ) 
dctM(hi, . . . ,h r ) 



(86) 
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Moreover, choosing the set of linearly independent vectors hi, . . . , h r so that 
h r = 27 is the gradient of the normalization constraint (Trp = 7 • h r /2), and 
defining the "mean" functionals (Hj) = TrpHj = hj ■ h r /A [in the paper, we 
assume r — 2, with Hi = H , the Hamiltonian operator, and Hi = / the identity, 
so that hi = H' = 2-/H and (Hi) = TrpH = 7 • H'/2 = (2 7 ) • (2 7 iT)/4 = 
hi • h 2 /4], (B) = b • h r /4 [in the paper, B = S, b = S'], and the "deviation" 
vectors Ahj = (hj — h r (Hj))/2, Ab = (b — h r (B))/2, it is easy to show that 

detM(b,hi,...,h») = detM(Ab,Ahi,...,Ah n _i) 
detM(hi,...,h„) ~ detM(Ahi,...,Ah„_i) ' 1 ' 

As a final remark, we write the following generalized form of the Cauchy- 
Schwarz inequality 

detM(Ab,Ahi,...,Ah„_i) < det M (Ahi, . . . , Ah„_i) Ab • Ab . (88) 

References 

[1] H.D. Zeh: The Physical Basis of the Direction of Time, Springer, 2001. 

[2] J. Uffink: in The Handbook of Philosophy of Physics, J. Butterfield and J. 
Earman, editors, Elsevier, 2006. 

[3] A. Bassi and G.C. Ghirardi: Physics Reports 379, 257 (2003). 

[4] D.Z. Albert: Brit. J. Phil. Set. 45, 669 (1994). D.Z. Albert, Time and 
Chance, Harvard University Press, 2001. 

[5] A. Kossakowski: Bull. Acad. Sci. Math. 20, 1021 (1972); A. Kossakowski: 
Rep. Math. Phys. 3, 247 (1972); R.S. Ingarden and A. Kossakowski: Ann. 
Phys. (N.Y.) 89, 451 (1975); G. Lindblad: Comm. Math. Phys. 48, 119 
(1976). 

[6] G.P. Beretta: Sc.D . thesis, MIT, Cambridge, MA, 1981; 
|arXiv:quant-ph/0509116| 

[7] V. Gorini, A. Kossakowski and E.C.G. Sudarshan: J. Math. Phys. 17, 821 
(1976). 

[8] W. Thomson, lord Kelvin: Proc. Royal Soc. Edinburgh 3, 139 (1852). 

[9] G.P. Beretta, E.P. Gyftopoulos, J.L. Park, and G.N. Hatsopoulos: Nuovo 
Cimento B 82, 169-191 (1984). 

[10] G.P. Beretta, E.P. Gyftopoulos, and J.L. Park: Nuovo Cimento B, 87, 
77-97 (1985). 

[11] G.P. Beretta: Found. Phys. 17, 365-381 (1987). 



31 



[12] G.P. Beretta: [arXiv:quant-ph/0612215| (2006); G.P. Beretta: Int. J. Quan- 
tum Information 5, 249 (2007). 

[13] J. Maddox: Nature 316, 11 (1985). 

[14] G.P. Beretta, A.F. Ghoniem, and G.N. Hatsopoulos, Editors: Meeting the 
Entropy Challenge, AIP Conf. Proc. 1033, 2008. 

[15] G.N. Hatsopoulos and E.P. Gyftopoulos: Found. Phys. 6, 15-31, 127-141, 
439-455, 561-570 (1976). 

[16] G.P. Beretta: Mod. Phys. Lett. A 20, 977 (2005). 

[17] G.P. Beretta: Phys. Rev. E 73, 026113 (2006). 

[18] G.P. Beretta: in Computer- Aided Engineering of Energy Systems, R.A. 
Gaggioli, Editor, (ASME Book H0341C, ASME, New York, 1986), pp. 
129-134. G.P. Beretta: in Second Law Analysis of Thermal Systems, M.J. 
Moran and E. Sciubba, Editors, (ASME Book 100236, ASME, New York, 
1987), pp. 17-24. G.P. Beretta: in Second Law Analysis of Heat Transfer in 
Energy Systems, R.F. Boehm and N. Lior, Editors, (ASME Book G00390, 
HTD 80, ASME, New York, 1987), pp. 31-38. 

[19] G.P. Beretta: Entropy 10, 160-182 (2008). 

[20] G.P. Beretta: Mod. Phys. Lett. A 21, 2799-2811 (2006). 

[21] W. Jones: J. Phys. C 15, 6597 (1982). 

[22] M. Courbage and I. Prigogine: Proc. Natl. Acad. Sci. USA, 80, 2412 (1983). 
[23] R.C. Dewar: J. Phys. A, 38 L371 (2005). 

[24] H. Struchtrup and W. Weiss: Phys. Rev. Lett., 80, 5048 (1998). 

[25] R.D. Levine and M. Tribus, Editors: The Maximum Entropy Formalism; 
The M.I.T. Press: Cambridge, Mass., 1979. 

[26] Y. Alhassid, N. Agmon and R.D. Levine: Phys. Rev. C 79, 1775-1788 
(1978). 

[27] R.D. Levine: J. Chem. Phys. 65, 3302-3315 (1976). 

[28] J.C. Keck: in The Maximum Entropy Formalism, R.D. Levine and M. 
Tribus, Editors; The M.I.T. Press: Cambridge, Mass., 1979. pp. 219-245. 

[29] G.P. Beretta and J.C. Keck: in Computer- Aided Engineering of Energy 
Systems, Second Law Analysis and Modeling, R.A. Gaggioli, Editor, ASME 
Book H0341C-AES, 1986, Vol. 3, pp. 135-139. 



32 



[30] G.P. Beretta: in Frontiers of Nonequilibrium Statistical Physics, Proc. 
NATO ASI, Santa Fe, 1984, G.T. Moore, G.T. and M.O. Scully, Editors, 
NATO ASI Series B: Physics 135, Plenum Press, New York, 1986, pp. 
193-204 and 205-212. 

[31] G.P. Beretta: in The Physics of Phase Space, Y.S. Kim and W.W. Zachary, 
Editors; Lecture Notes in Physics, Vol. 278, Springer- Verlag, pp. 441-443 
(1986). 

[32] H. Ozawa, A. Ohmura, R.D. Lorenz, and T. Pujol: Reviews of Geophysics, 
41, 1018 (2003). 

[33] H.C. Ottinger and M. Grmela: Phys. Rev. E 56, 6620, 6633 (1997). 

[34] A. Caticha: in Bayesian Inference and Maximum Entropy Methods in Sci- 
ence and Engineering: 21st Intl. Workshop, R. L. Fry, Editor; AIP Conf. 
Proc. 617, 302 (2002). 

[35] B.J. Chung and A. Vaidya: Physica D 237, 29452951 (2008). 

[36] L.M. Martyusheva and V.D. Seleznev: Physics Reports 426, 1-45 (2006). 

[37] S. Gheorghiu-Svirschevski: Phys. Rev. A 63, 022105 (2001); ibid, 63, 
054102 (2001). 

[38] M. Lemanska and Z. Jaeger: Physica D 170, 72 (2002). 

[39] A. Caticha: in Bayesian Inference and Maximum Entropy Methods in Sci- 
ence and Engineering: 20th Intl. Workshop, A. Mohammad-Djafari, Editor; 
AIP Conf. Proc. 568, 72 (2001). 

[40] E.P. Gyftopoulos and G.P. Beretta: Thermodynamics. Foundations and 
Applications, Dover, 2005 (first edition, Macmillan, 1991). 

[41] W.K. Wootters: Phys. Rev. D 23, 357-362 (1981). 

[42] P. Salamon, J.D. Nulton, and R.S. Berry: J. Chem. Phys. 82, 2433-2436 
(1985). 

[43] S.L. Braunstein and CM. Caves: Phys. Rev. Lett. 72, 3439-3443 (1994). 
[44] G.P. Beretta: [arXiv:quant-ph/0112046| (2001). 

[45] A. Messiah: Quantum Mechanics, Wiley, 1976, p. 320; L. Mandelstam and 
I.G. Tamm, J. Phys. USSR 9, 249 (1945). 

[46] G.P. Beretta: [arXiv:quant-ph/0511091| (2005). 

[47] G.P. Beretta: Int. J. Theor. Phys. 24, 119 (1985). 

[48] G.N. Hatsopoulos and J.H. Keenan: Principles of General Thermodynam- 
ics, Wiley, 1965. 



33 



[49] M. Ozawa: Annals of Physics 311, 350 (2004). 

[50] J.L. Park and H. Margenau: Int. J. Theor. Phys. 1, 211 (1968). 

[51] R. Alicki, M. Horodecki, P. Horodecki, and R. Horodecki: Open Systems 
and Information Dynamics 11, 205 (2004). 

[52] W. Band and J.L. Park: Found. Phys. 1, 133 (1970), 1, 211, 339 (1971); 
U. Fano: Rev. Mod. Phys. 29, 74 (1957). 

[53] A.E. Allahverdyan, R. Balian, and Th.M. Nieuwenhuizen: Europhysics Let- 
ters 67, 565 (2004). 



34 



